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Abstract. We use a time-dependent dynamical mean-field-hydrodynamic model to predict and study bright 
solitons in a degenerate fermion-fermion mixture in a quasi-one-dimensional cigar-shaped geometry using 
variational and numerical methods. Due to a strong Pauli-blocking repulsion among identical spin-polarized 
fermions at short distances there cannot be bright solitons for repulsive interspecies fermion-fermion inter- 
actions. However, stable bright solitons can be formed for a sufficiently attractive interspecies interaction. 
We perform a numerical stability analysis of these solitons and also demonstrate the formation of soliton 
trains. These fermionic solitons can be formed and studied in laboratory with present technology. 

PACS. 03.75.Ss Degenerate Fermi gases - 05.45.Yv Solitons 



Recent observations |l|2|3|4j and associated experimen- 
tal [5I5I7j and theoretical |8I9I10I11I12] studies of a degen- 
erate Fermi gas (DFG) by sympathetic cooling in the pres- 
ence of a second boson or fermion component suggest the 
possibility of soliton formation. Apart from the observa- 
tion of a DFG in the following degenerate boson-fermion 
mixtures (DBFM) '^•^Li [3], ^^Na-^Li [4j and *^^Rb-^"K 
|5l6j , there have been studies of degenerate spin-polarized 
fermion-fermion mixtures (DFFM) "^"K-^^K [J and ^Li- 

6Liia. 

Bright solitons in a Bose-Einstein condensate (BEG) 
are formed due to an attractive nonlinear atomic inter- 
action [13]. As the interaction in a pure DFG at short 
• distances is repulsive due to strong Pauli blocking, there 
cannot be bright solitons in a DFG. However, it has been 
, demonstrated ^14)15j that bright solitons can be formed 
■ in a DBFM in the presence of a sufficiently strong boson- 
fermion attraction which can overcome the Pauli repulsion 
among identical fermions. 

We demonstrate the formation of stable fermionic bright 
, solitons in a DFFM for a sufficiently attractive interspecies 
' fermion-fermion interaction. In a DFFM, the coupled sys- 
tem can lower its energy by forming high density regions, 
the bright solitons, when the attraction between the two 
types of fermions is large enough to overcome the Pauli re- 
pulsion among identical fermions. We use a coupled time- 
dependent mean-field-hydrodynamic model for a DFFM 
and consider the formation of axially-free localized bright 
solitons in a quasi-one-dimensional cigar-shaped geome- 
try using numerical and variational solutions. The present 
model is inspired by the success of a similar model sug- 
gested recently by the present author in the investigation 
of collapse [12] and bright [15] and dark [16] solitons in 
a DBFM. We study the condition of modulational in- 



stability of a constant-amplitude solution in this model 
and demonstrate the possibility of the formation of bright 
solitons. We also present a numerical stability analysis of 
these robust bright solitons and consider the formation of 
a soliton train in a DFFM by a large sudden jump in the 
interspecies fermion-fermion scattering length near a Fes- 
hbach resonance, experimentally observed in both 
and 40K-40K [17]. 

We use a simplified mean-field-hydrodynamic Lagran- 
gian for a DFG used successfully to study a DBFM |12I15I16| . 
The virtue of the mean-field model over a microscopic de- 
scription is its simplicity and predictive power. To develop 
a set of time-dependent mean-field-hydrodynamic equa- 
tions for the interacting DFFM, we use the following La- 
grangian density [12|15j 
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where j — 1,2 represents the two components, Tpj the 
complex probability amplitude, rij — IV'jP the real proba- 
bility density, * denotes complex conjugate, mj the mass, 
Aj — ?i^(67r^)^/^/(2mj), the interspecies coupling gi2 — 
2ttTi ai2/mR with rUR = mim2/(mi -I- TO2) the reduced 
mass, and ai2 the interspecies fermion-fermion scattering 
length. The number of fermionic atoms Nj is given by 
J dxrij (r) — Nj . The trap potential with axial symmetry 
is Vj{r) — ^mjUj'^{p'^ + v^z^') where lo and vuj are the 
angular frequencies in the radial (p) and axial (z) direc- 
tions with V the anisotropy. The interaction between iden- 
tical intra-species fermions in spin-polarized state is highly 

5 /3 

suppressed due to Pauli blocking terms SAjUj /5 and 
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has been neglected in Eq. ([T]). The kinetic energy terms 
h I VrV'j P/CSTij) in Eq. ([T]) contribute httle to this prob- 
lem compared to the dominating Pauli-blocking terms. 
However, its inclusion leads to an analytic solution for 
the probability density everywhere fl5J . 

With the Lagrangian density |T]), the following Eulcr- 
Lagrange equations can be derived in a straight-forward 
fashion pITS] : 



ih 



d 



dt 6m i 



A 2/3 



gunk 



= 0, (2) 
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where j ^ k — 1,2. This is essentially a time-dependent 
version of a similar time-independent model used recently 
For large nj, both lead to the Thomas- 



[(Mj - Vj)/^3?^^ [I0TT2] with Hj the 



for fermions [10] 
Fermi result nj 
chemical potential. As the bright solitons of this rapid note 
are stationary states, they could be obtained by the time- 
independent approach used in Ref. [10]. The stationary 
approach of Ref. [10] has passed rigorous tests of compar- 
ison of the hydrodynamic-mean-field spectra of localized 
fermions with the spectra calculated in the coUisionlcss 
regime within the random-phase approximation (RPA). 
The results of mixing-demixing and collapse of the hydro- 
dynamic approach arc in agreement with the RPA anal- 
ysis [TT]. The detailed behavior of collective excitation of 
trapped fermions has also been found to agree with that 
obtained by an RPA analysis [TU]. For a description of 
stationary solitons (e.g., of Fig. 1) we could have used the 
well-established formulation of Ref. JSH to obtain identi- 
cal results, as the present time-dependent dynamical de- 
scription and the time-independent approach of [10] yield 
identical results for stationary states. However, we shall be 
using the present time-dependent dynamical formulation 
to study the nonequilibrium generation soliton trains, in 
addition. 

We reduce three-dimensional Eqs. ([2]) to a minimal 
quasi-one-dimensional form in a cigar-shaped geometry 
with V << 1, where the radial motion is frozen in the 
ground state of the harmonic trap and the dynamics is car- 
ried by the axial motion. For radially-bound and axially- 
free solitons we eventually set i' — 0. Following Ref. 15] 
this reduction can be done in a straight-forward fashion 
and we quote the final results here: 



,4/3 



^,(y,r) = 0, (3) 



1,2 represents the two solitons, r — 
= 9(6^iVj)2/3/5, Njk = 12\ai2\Nk/l, 
3toi — 3m2. Here we employ 



d 

where ^ k = 
tuj/2, y = z/l, Njj 
I = ^/h/{u)m), with m 
equal-mass fermions, a negative ai2 corresponding to at- 
traction, and normalization \4>i{y-: 'T)\'^dy = 1- In Eqs. 
((3|) a sufficiently strong attractive fermion-fermion cou- 
pling A^jfc|<^feP(j 7^ k) can overcome the Pauli repulsion 
Njj\(f)j\'^/^ and form bright solitons. 

Now we perform a stability analysis of constant-ampli- 
tude solutions of Eqs. ((3|) and study the possibility of gen- 
eration of solitons in the symmetric case: iVi — N2, when 



Fig. 1. (Color online) The solitons |9ij(y)| of Eq. ([3]) vs. y (in 

dimensionless units) for A^i = 44, N2 = 56, ai2 = —0.3 nm, 
while TVii 159, Nx2 ~ 203, N2x ^ 160, and N22 ~ 187. 
The variational (ifv) and numerical (ip) solutions of Eq. Q for 
A^i = = 50 are also shown. 
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Fig. 2. (Color online) The propagation of fermionic solitons 
(a) \(l)i{z,t)\ and (b) \4>2{z,t)\ of Fig. 1 vs. z and t. At t = 100 
ms (marked by arrows) the bright solitons are set into small 
breathing oscillation by suddenly changing Ni = 44 and = 
56 to N2 = 44, iVi = 56. 



bi ~ (f)2 = (p and these equations reduce to 
d 



4/3 , I |2 



(p(2/,r) = 0, (4) 



= N21- We consider 
= Aoexp(i(5) = 

PAy^r)] of Eq. Q under small pertur- 



where {3 = Nn — N22 and 7 = N12 
the constant-amplitude solution [1^ 

Aq exp[i(7Agr 

bation: = {Aq + A) exp(«i5), where A = A{y, r) and Aq 
the amplitude. Substituting this perturbed solution in Eq. 
(U, and for small perturbations retaining only the linear 
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terms in A we get 



OA d^A 



dr dy^ 



-pAt^'iA + A*) + 7A2(A + A*) = 0. (5) 



We consider the plane-wave perturbation A{y, r) — Ai- 
cos{Kt — fly) + iA2 sm{KT — fly) in Eq. ([5]). Then sep- 
arating the real and imaginary terms and eliminating Ai 
and A2 we obtain the dispersion relation K — ±i7[J7^ 

-{^-yAl - 4/3Ao/^/3)]^/2. For stability of the plane- wave 
perturbation, K has to be real. This happens for 27^^ < 

AI3AI'^/?, or -iAI'^ < 2/3/3. However, K can become imag- 
inary for 7^9^"^ > 2/3/3 and the plane-wave perturbations 
can grow exponentially with time t. This is the domain of 
modulational instability of a constant-intensity solution, 
signalling a tendency of spatially localized bright solitons 
to appear. We also performed this analysis in the case of 
non-symmetric coupled equations ([3]) and quote the result 

2/3 2/3 

here. The condition for instability is iVi2A'2i^iQ > 
4AfiiA''22/9 [H], where ^10 and A20 are the amplitudes of 
the two solutions. 

Next we present a variational analysis of Eq. (|4]) based 
on the normalized Gaussian trial wave function [20j fv{y, t) 
= v/l/[a(T)0F] exp[-2/V{2a2(r)} + ih{T)y^ /2], where a 
is the width and h the chirp. The Lagrangian density for 
Eq. (HI is the one-term version of Eq. (jT]) , which is evalu- 
ated with this trial function and the effective Lagrangian 
L = C{ify)dy becomes 



L = 



b + 



2b^ - 



V2 7 12V3 (3 



/^a3 5^5 ^i/3a8/3 



(6) 



The variational Euler-Lagrangian equations for a and b 
can then be written and solved in a standard fashion [5D] 
to yield the differential equation for the width: (Pa/dr^ — 
[4-a7A/27^-Fa4/3(8/3\/3)/(57ri/375)]/a3^ .pj^g variational 

result for width a follows by setting the right hand side 
of this equation to zero, from which the variational profile 
for the soliton can be obtained [20J. 

We solve Eqs. ^ for bright solitons numerically us- 
ing a time-iteration method based on the Crank-Nicholson 
discretization scheme elaborated in Ref. [21] using time 
step 0.0002 and space step 0.015. We perform a time evolu- 
tion of Eqs. ([3]) introducing an harmonic oscillator poten- 
tial y^ and setting the nonlinear terms to zero, and start- 
ing with the eigenfunction of the linear harmonic oscillator 
problem. The extra harmonic oscillator potential, which 
will be set equal to zero in the end, only aids in starting 
the time evolution with an exact analytic form. During the 
time evolution the nonlinear terms are switched on and the 
harmonic oscillator potential is switched off slowly and the 
time evolution continued to obtain the final converged so- 
lutions. 

In our numerical study we take I = 1 /im and consider 
a DFFM consisting of two electronic states of ""^K atoms. 
This corresponds to a radial trap of frequency w w 27r x 83 
Hz. Consequently, the unit of time is 2/uj k A ms. 

First we solve coupled Eqs. with Ni = 44, N2 = 56, 
and ai2 = —0.3 nm. The soliton profile in this case is 
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Fig. 3. (Color online) The propagation of fermionic solitons (a) 
|(/>i(z,t)| and (b) |</>2(z,t)| vs. z and t. At t = 100 ms (marked 
by arrows) the same initial solitons of Figs. 1 are set into small 
breathing oscillation by suddenly changing 0i l.l(/>i and 

4>2 — » 0.9(/)2. 



shown in Fig. 1, where we also plot the variational and 
numerical solutions of Eq. @ for Ni = N2 — In this 
case in Eq. @ (3 ~ 173 and 7 ~ 181.6 leading to a vari- 
ational width a « 0.043 and a variational soliton profile 
(pi, « 3.62 exp(— 270?/^). The variational result agrees well 
with the numerical solutions. 

At this stage it is pertinent to see if Friedel oscillations 
[22] of density of the localized fermions are small so that 
the effective description of the solitons is valid. An one- 
dimensional degenerate Fermi gas of N atoms filled up 
to Fermi sea has a spatial extension 2Lp = l\/2N — 1, 
where I is a measure of confinement [23'. In the pres- 
ence of an harmonic trap, I is the harmonic trap length 
and is smaller than the spatial extension of the confined 
fermions. In the case of the soliton of Fig. 2, a typical mea- 
sure of / could be 0.05 /^m so that for about 50 fermions 
considered here 2Lp = 0.5 /im. The Fermi momentum 
of an one-dimensional Fermi gas \s kp = ttN/{2Lp) [23j . 
The spatial wave-length of Friedel oscillation is [23] A = 
Tr/kp = 2Lp/N « 0.01 /im, much smaller than the soli- 
ton width of 0.1 /im. This qualitative analysis resulting in 
small Friedel oscillation supports the effective description 
used this rapid note for a DFFM. 

To test the robustness of these solitons we infiicted 
different perturbations on them and studied the resultant 
dynamics numerically. First, after the formation of the 
solitons we suddenly changed the fermion numbers from 
Ni = 44 and N2 = 56 to Ni = 56, N2 = 44 at time t = 100 
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Fig. 4. (Color online) Soliton trains of two and three solitons 
formed upon removing the harmonic trap and jumping the 
nonlinearities at t = from Njj ~ 173, Njk ~ 145, fc 7^ j = 1, 2 
to (a) Njj « 173, Njk « 273, and to (b) iV^^ ^ 173, Njk ^ 309,, 
respectively. 

ms. This corresponds to a sudden change of nonUnearities 
from iVii w 159, N12 ~ 203, N21 « 160, and N22 « 187 
to iVii « 187, N12 « 160, N21 « 203, and N22 ^ 159. The 
resultant dynamics is shown in Figs. 2 (a) and (b). Due to 
the sudden change in nonlinearities the fermionic bright 
solitons are set into stable non-periodic small-amplitude 
breathing oscillation. Next on the same initial solitons of 
Figs. 2, at i = 100 ms, we suddenly inflict the perturba- 
tion 01 1.101 and (f>2 0.902 and follow numerically 
the time evolution. The result of simulation is shown in 
Figs. 3 (a) and (b). We find that the solitons again con- 
tinue non-periodic breathing oscillation and stabilize at 
large times. We also gave a small displacement between 
the centers of these solitons. We find that after oscillation 
and dissipation the solitons again come back to the stable 
configuration. 

During the time evolution of Eqs. ^ if the nonlinear- 
ities are changed by a small amount or changed slowly, 
usually one gets a single stable soliton when the final non- 
linearities are appropriate. However, if the nonlinearities 
are jumped suddenly by a large amount, a soliton train 
can be obtained as in the experiment with BEC jT3]. To 
illustrate this we consider the solution of Eqs. ^ with 
nonlinearities Njj « 173 and Njk « 145, j ^ fc = 1, 2 and 
harmonic oscillator trap y^. After the formation of the 
solitons we suddenly jump the off-diagonal nonlinearities 
to Njk ~ 273 and also switch off the harmonic trap at 
time t = 0. Then after some initial noise and dissipation 



the time evolution of Eqs. ([3|) generates two slowly reced- 
ing bright solitons of each component as shown in Fig. 4 
(a). More solitons can be generated when the jump in the 
nonlinearities is larger. In Fig. 4 (b) we show the gener- 
ation of three slowly receding solitons of each component 
upon a sudden jump of the off-diagonal nonlinearities to 
Njk ~ 309 from the same initial state as in Fig. 4 (a). The 
formation of soliton trains from a stable initial state is due 
to modulational instability 18J. The sudden jump in the 
off-diagonal nonlinearities could be effected by a jump in 
the interspecies scattering length ai2 obtained by manipu- 
lating a background magnetic field near a fermion-fermion 
Feshbach resonance [17] . 

In conclusion, we use a coupled mean-field-hydrodyna- 
mic model for a DFFM to study the formation of bright 
solitons and soliton trains in a quasi-one-dimensional ge- 
ometry by numerical and variational methods. We find 
that an attractive interspecies interaction can overcome 
the Pauli blocking repulsion and form fermionic bright 
solitons in a DFFM. The stability of the present solitons is 
demonstrated numerically through their sustained breath- 
ing oscillation initiated by a sudden small perturbation. 
We also illustrate the creation of soliton trains upon a 
sudden large jump in off-diagonal nonlinearities. Bright 
solitons and soliton trains have been created experimen- 
tally in attractive BECs in the presence of a radial trap 
only without any axial trap fl3]. In view of this, fermionic 
bright solitons and trains can be created in laboratory in 
a DFFM in a quasi-one-dimensional configuration. Here 
we used a set of mean-field equations for the DFFM. A 
proper treatment of a DFG or DFFM should be done using 
a fully antisymmetrized many-body Slater determinant 
wave function |8ll4j as in the case of scattering involv- 
ing many electrons |24j . However, in view of the success of 
a fermionic mean-field-hydrodynamic model in studies of 
collapse 12 , bright [15] and dark solitons [16] in a DBFM, 
and of mixing-demixing [25j and black solitons [26 in a 
DFFM, we do not believe that the present study on bright 
solitons in a DFFM to be so peculiar as to have no general 
validity. 

The work is supported in part by the CNPq and FAFESF 
of Brazil. 
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